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Abstract 



The pair- and three-proton interaction potentials for metal-phase hydrogen 
are calculated. Irreducible three-proton interactions are shown to be essential 
in the development of the metal hydrogen structure. Possible manifestations 
of the three-ion interactions in the structure of quickly quenched metallic 
systems are discussed. Existence of amorphous metals is related to mani- 
festations of the three-ion interactions in non-equilibrium conditions when 
amorphisation front travels through a liquid metal. 

I. INTRODUCTION 

The inclusion of many-particle interactions into consideration of equilibrium and kinetics 
properties is the urgent problem in the theory of condensed matter PL O. Calculation of 
the many-particle interaction potentials is the key to this problem. For simple dielectric 
and metal liquids and for crystals pair interaction potentials are known at least at a phe- 
nomenological level, about the three-particle interaction it cannot be said. This situation is 
explained by general lack of development of the theory of highly nonuniform electron gas. 
Thus to calculate the three-particle interactions in various condensed systems one needs to 
know in an explicit form the three-point nonlinear response functions involving cooperative 
phenomena in the many-electron system, which is formed from atomic electron shells. Nowa- 
days the problems of calculation of such functions for the majority of systems are not even 
formulated. It should be stressed that for the most systems even the polarization operator 
of nonuniform electron gas, which makes up these systems, remains unknown. But there 
are same exceptions. This refers to the simple metals which can be closely approximated 
by the uniform gas of free electrons. A discrete nature of the ion subsystem can be taken 
into account by perturbation theory. The nonlinear response function of uniform gas cor- 
responding to three- and four-pole diagrams have been calculated within the framework of 
many-particle theory of metals [[|. 

Calculations of the structural and thermodynamical equilibrium volume properties of 
metal-phase hydrogen point up the significance of inclusion of many-particle ion interac- 
tions, and suggest the possibility of a metastable metallic state at zero external pressure 
0. Comparison of energies corresponding to different crystal structures have led to the 
conclusion of liquidity tendencies in the structure of metal hydrogen. Energy difference of 
crystal and liquid phases approaches zero with decreasing pressure ||. So the liquid state 
of metal hydrogen is possible to form at zero external pressure. 

Calculations of energy and other properties of crystal metals are carried out in the 
reciprocal space and do not require the explicit form of interaction potentials of ion clusters. 
In amorphous, liquid and heterogeneous metals the calculations in configurational space are 
preferable. It provides insight into the nature of local atomic order and gives the assessment 
of correlations in the positions of atom clusters. 

The occurrence of a maximum in the temperature dependence of the third virial coeffi- 
cient evidences directly the existence of three-particle interactions in simple dielectric liquids 
IHII . Taking into account the three-atom interaction in the asymptotic form of the Akselrod- 
Teller three-dipole interaction potential requires the inclusion of a cutoff parameter (the 
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characteristic length) into the theory and does not explain the wave-vector dependence of 
the third virial coefficient (on the example of krypton) || . 

In metals the many-ion interactions are directly proven to exist by the anomalies in 
phonon spectra which could not be reduced to the Kohn's ones, and conform with the 
singularities of the many-pole diagrams of uniform electron gas [Q]. 

Interaction potential in liquid Na and K have been calculated previously || for three- 
ion equiangular configurations. The pair-interaction potentials for ions in metal hydrogen 
have been calculated in j7|, where the third-order contributions with respect to electron-ion 
interaction were shown to be of fundamental importance for development of attractive part 
of potential. Pair interactions in a liquid metal hydrogen have been calculated within the 
context of the density-functional method |§. The potentials were formulated in terms of 
direct correlation functions of electron-ion system. 

Here we present the pair- and three-ion interaction potentials in metal hydrogen at the 
Wigner-Seitz radius r s = 1, 65 calculated on a basis of the many-particle theory of metals 
within the third-order perturbation theory. The computed array of values of three-ion 
interaction potential makes possible the further calculations of the contributions of these 
interactions into various properties of metal hydrogen. Possible outcomes of the three-ion 
interactions affecting local atomic order in amorphous metal systems are discussed. 



II. CALCULATION RESULTS 

In the context of many-particle theory of non-transition metals the energy of a metal is 
calculated by the use of the adiabatic approximation for electron-ion system. The electron 
gas energy E e , provided that ions positions are fixed, can be calculated with the theory of 
perturbations in the potentials of electron-ion and electron-electron interactions ||. The 
energy E e fulfills the role of an electron term upon the treatment of the properties of ion 
subsystem in metal. Together with kinetic energy of ions and energy of their direct interac- 
tion, the energy E e comprises the effective Hamiltonian of the ion subsystem of metal. E e 
can be considered as a sum of contributions independent of ion positions, dependent on the 
locations of separate ions, ion pairs, triplets, etc. [[J: 

E e = V?0 + X)Vl( R w) + ^|E^(Rn,R m ) + ^ E V9 3 (R n ,R m ,R;) H (1) 

Each term of the series (1) describes interactions of ion groups through the surround- 
ing electron gas and can be represented as a power series in the potential of electron-ion 
interaction ||: 

oo 

^(R 1; R 2 ) = ^$f ) (Ri,R 2 ), 

i=2 
oo 

^ 3 (Ri,R2,R 3 ) = £$?(Ri,R2,R 3 ), (2) 

i=3 

etc., where $^(Ri, . . . , R n ) represents the indirect interaction of the n ions through electron 
gas in the /c-order perturbation theory in electron-ion interaction. 



3 



The indirect interaction of two ions described by the potential $ 2 (|Ri — R2I) is well 
studied. The potential $2 (|R-i — R2I), corresponding to the energy of indirect two-ion 
interaction in the third-order in potential of electron-ion interaction, has been calculated 
previously in |J for liquid sodium and potassium, and also for metal hydrogen JF]. It is not 
difficult to obtain an expression for the $ 2 3) (|Rx - R 2 |): 

00 00 1 

^2 } (R) = 4^4 / dq X ql J dq 2 q\ J dxV( gi ) V(q 2 )V(q 3 ) x 
00-1 



v r (3)/ ^m( qi R) 

x r w (gi,g 2 ,g 3 ) R , (3) 

where V(g) is the form-factor of the electron-ion interaction potential; r^ 3 ^(gi, q 2 , q 3 ) is the 
sum of three-pole diagrams. 

In the third-order perturbation theory in V(q) the indirect three-ion interaction (the 
irreducible three-ion interaction) is defined by the expression 

00 00 
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$S 3) (i*i, R 2 , R 3 ) = —J dq x ql J dq 2 q 2 2 x 



1 

x / dzV(q 1 )V(q 2 )V(q 3 )T^(q 1 ,q 2 ,q 3 ) J dxx 





x cos x[ q±Ri 



p2 j_ p2 p2 

2R\R 2 



+ g 2 ^2^ x 



X 



x Jo ^2 (l - x 2 ) V2 - z 2 ^j j , (4) 

where Jq{x) is the Bessel function of zero order; z = cos(qi, q 2 ); (fe = (qf + q 2 + 2qiq 2 x) 1 ^ 2 ; 
Ri, R 2 , R 3 - the distances between the vertices of a triangle formed by the protons. 

The pair- and three-proton interaction potentials were calculated at the Wigner-Seitz 
radius r$ = 1, 65 which corresponds to zero pressure in the zeroth model of a metal. A 
permittivity function in the Heldart-Vosko form was employed. 

Figure 1 shows computed potentials of the two-proton interactions e 2 /R + ^^(R), 
<S> { 2 \R), and <p*(R) = e 2 /R + $ { 2 \R) + $ 2 \R). It is obvious that the interaction 
significantly renormalizes potential e 2 /R + $ 2 2 (R)- Noteworthy is a minimum in the repul- 
sive part of the potential <p*(R), which arises as a result of the contribution $2 (i?). At 
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r$ = 1, 72 this minimum becomes deeper and as a result the potential <f*{R) takes the form 
typical for simple metals. At r$ < 1,65 the minimum in the repulsive part of <p*{R) turns 
shallow and its position shifts towards smaller R. 

The potentials of irreducible three-proton interaction $3 (Ri, R%, R3) calculated at dif- 
ferent lengths of sides of a three-proton triangle are listed in Table. It should be noted 
that the number of various (R\, R2, R3) sets is restricted by the triangle axiom to which the 
distances Ri,R2 and R3 are subject. 

The calculation results can be presented more conspicuously and informatively as plots 
of versus geometry parameters for selected configurations of protons. Thus, Figure 2 demon- 
strates the potential curve <f>f\R,R,R), that is, one for proton equiangular triplets. The 
potential $3 (R, R, R) corresponds to strong attraction of the triplet at short distances and 
oscillates at the long ones. 

Figure 3 shows plots of the potential <3> 3 3 ' ) for isosceles ion triplets. It is clear that 
attractability of proton triplet is actually short-range. An equilibrium local atomic order 
results from the balance of direct proton interactions, the indirect two-proton ones, and of 
the irreducible three-proton interactions. 

Figure 4 presents a proton energy plotted against the distance to the two others (the 
isosceles proton triplet). Two protons are placed in abscissa axis at the points (1,0,0) and 
(-1,0,0) (distances in atomic units). The third proton is sited in ordinate axis. In Fig. 4 

(3) 

is shown that inclusion of the irreducible potential $ 3 is of fundamental importance for 
determination of the equilibrium local order in proton spacing. Attractive nature of the $ 3 
potential have to reduce the average interionic distances in equilibrium state of metal. 



III. DISCUSSION AND CONCLUSION 

A comprehensive analysis of the atomic properties of simple metals leads to the conclusion 
that of basic importance is the consideration of three-ion interactions as in polyvalent metals 
||, and also in metal hydrogen where the electron- ion interaction is free from a non-coulomb 
part ||. Calculations of the two-ion interaction potentials for various simple metals display 
the similarity of their behavior in all such metals ||10|| . Our calculations $ 3 3 ' ) for equilateral 
proton triplets, along with the results of |J , demonstrate the similarity of the $ 3 3 ^ behavior 
in simple metals. 

Results of our calculations suggest certain inferences about the possible manifesta- 
tions of three-particle interactions through the structure of metal systems. The potential 
$ 3 (Ri, R%, R3) is symmetric about interchange of ions and therefore in equilibrium condi- 
tions the most probable three-ion configurations in liquid metals would be the equilateral 
triplets. 

Experiments on the quick quenching of metals and alloys show that amorphous phase is 
readily obtainable in the polyvalent metals and their alloys [|TTJ where consideration of many- 



ion interactions are essential. Under the quick quenching conditions aluminum alloys form 
quasicrystalline structures with specific 5-fold axes of symmetry forbidden for crystals []T2 



The possible reason of quasicrystal formation is manifestation of the three-ion interactions 
as the amorphisation front moves through a liquid metal. 

When two of the three ions are fixed on a surface or located within amorphous phase and 



5 



the third ion remains in the liquid, energetically advantageous would be an isosceles triplet 

(3) 

rather than equilateral one. Owing to the attractive nature of $3 the third ion shifts so 
that one of the angles in the triplet will be more than 60°. Such isosceles triplets can serve 
as a basis for generation of amorphons in an amorphous phase. Existence of the 5-fold axes 
of symmetry in quasicrystals corresponds to the presence of the isosceles triplets with an 
angle of 72°. Thus one of the manifestation of three-ion interactions can be the creation of 
amorphons while propagating the amorphization front in liquid metal. 

Acknowledgement 1 The authors are grateful to Prof. Yu. P. Krasny for stimulating 
discussion. 
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0.00029 
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-0.00028 
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0.00065 
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0.00038 
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6 


-0.00024 
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-0.00033 
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-0.00134 
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-0.00030 
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0.00048 
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0.00006 
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0.00027 
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4 
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0.00025 
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-0.00019 
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-0.00007 
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-0.00021 
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0.00006 
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-0.00004 
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0.00007 
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0.00006 
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0.00003 
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(distances in the Table in atomic units, energy in Ry) 

Figure Captions 

Fig. I. Potential of the two-proton interactions: 1 - e 2 /R + $£ 2) (.R); 2 - ^(R); 
3 - <p*(R) = e 2 /R + $?\r) + $ { ?\R) 
Fig. 2. Potential $J 3) (R, R, R) 

Fig. 3. Potential <^ 3) : 1 - $f } (1 a.u., R, R); 2 - $jj 3) (2 a.u., R, R); 3 - $jj 3) (3 a.u., R, R) 
Fig 4. Potential ip(R) = 2<p*(R) + $^ 3) (2 a.u., R, R): 1 - 2<p*(R); 2 - $f } (2 a.«., R, R); 
3 - (see explanations in the text) 
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